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ABSTRACT 

We have used the Australia Telescope Compact Array (ATCA) at 95 GHz to carry out 
continuum observations of 130 extragalactic radio sources selected from the Australia 
Telescope 20 GHz (AT20G) survey. We use a triple-correlation method to measure 
simultaneous 20 and 95 GHz flux densities for these objects, and over 90% of our 
target sources are detected at 95 GHz. We demonstrate that the ATCA can robustly 
measure 95 GHz flux densities with an accuracy of ~10% in a few minutes for sources 
stronger than about 50 mJy. 

We measure the distribution of radio spectral indices in a flux-limited sample 
of extragalactic sources, and show that the median 20-95 GHz spectral index does 
not vary significantly with flux density for S20 > 150 mJy. This finding allows us to 
estimate the extragalactic radio source counts at 95 GHz by combining our observed 
20-95 GHz spectral-index distribution with the accurate 20 GHz source counts mea- 
sured in the AT20G survey. 

Our derived 95 GHz source counts at flux densities above 80 mJy are significantly 
lower than those found by several previous studies. The main reason is that most radio 
sources with flat or rising spectra in the frequency range 5-20 GHz show a spectral 
turnover between 20 and 95 GHz. As a result, there are fewer 95 GHz sources (by al- 
most a factor of two at 0.1 Jy) than would be predicted on the basis of extrapolation 
from the source populations seen in lower-frequency surveys. We also derive the pre- 
dicted confusion noise in CMB surveys at 95 GHz and find a value 20-30% lower than 
previous estimates. 

The 95 GHz source population at the flux levels probed by this study is domi- 
nated by QSOs with a median redshift z ~ 1. We find a correlation between optical 
magnitude and 95 GHz flux density which suggests that many of the brightest 95 GHz 
sources may be relativistically beamed, with both the optical and millimetre contin- 
uum significantly brightened by Doppler boosting. 

Key words: radio continuum: galaxies — galaxies: active — quasars: general — radio 
continuum: general — cosmic microwave background — techniques: interferometric 



1 INTRODUCTION 

The radio sky at 90-100 GHz remains largely uncharted ter- 
ritory. At frequencies above 5 GHz, traditional surveys of 
even modest areas of sky become impractical because of the 
small field of view of current telescopes. For example, the ef- 
fective field of view of the NRAO VLA telescope is 0.22 deg 2 
at 1.4 GHz, but less than 0.0002 deg 2 at its highest frequency 
of 50 GHz. Most radio telescopes designed to operate at mil- 



limetre wavelengths have similarly small fields of view* , so 
that blind surveys of large areas of sky are only possible by 
using specially-designed instruments or techniques. 

Figure 1 summarizes the sensitivity limits of large-area 
radio continuum surveys which cover at least a quarter of 

* In general, the area of sky covered by the primary beam of a 
radio telescope scales as v~ 2 for increasing frequency v. 
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Figure 1. Detection limits of the most sensitive large— area radio 
surveys available at frequencies from 74 MHz to 94 GHz. Only sur- 
veys which cover at least 25% of the sky have been included. The 
surveys plotted are the VLA Low-Frequency Sky Survey (VLSS, 
Cohen ct al. 2007), Wcstcrbork Northern Sky Survey (WENSS, 
Rcngclink et al. 1997), Sydney University Molonglo Sky Survey 
(SUMSS, Bock et al. 1999),NRAO VLA Sky Survey (NVSS, Con- 
don et al. 1998), single-dish PMN and GB6 surveys (Griffith & 
Wright 1993; Gregory et al. 1996), Australia Telescope 20 GHz 
Survey (AT20G, Ricci ct al. 2004) and the WMAP extragalactic 
source catalogue (Hinshaw et al. 2007). 



the sky. There are currently only two large-area surveys 
above 5 GHz. The WMAP Point Source Catalogue (Bennett 
et al. 2003; Hinshaw et al. 2007) detects sources down to 0.8- 
1.0 Jy in five frequency bands from 23-94 GHz. The AT20G 
survey (Ricci et al. 2004; Sadler et al. 2006) uses a wide- 
band (8 GHz bandwidth) analogue correlator to scan the sky 
rapidly at 20 GHz, detecting sources down to a limit of ~ 
40 mjy, and currently represents the highest radio frequency 
at which a blind survey of a large fraction of the sky has been 
carried to a detection limit below 100 mjy. 

There are several reasons why we would like to learn 
more about the extragalactic radio-source population at 
frequencies well above those probed by existing large-area 
surveys. Corrections for extragalactic foreground confu- 
sion (e.g. De Zotti et al. 2005) will be critical for next- 
generation surveys of the Cosmic Microwave Background 
(CMB) like the Planck mission, which aims to measure 
CMB anisotropies on angular scales of 5-30 arcmin in the 
frequency range 30-350 GHz. Such corrections require an ac- 
curate knowledge of the foreground point-source population 
at frequencies of 30 GHz and above. 

Observations at 95 GHz are particularly relevant to 
Planck, since the Planck 100 GHz channel and the nearby 
143 GHz channel are the "cleanest" ones for CMB studies in 
both temperature and polarization on small angular scales, 
where the main foreground contamination comes from extra- 
galactic sources. The 100 and 143 GHz channels will also be 
the most sensitive Planck channels^ (Lamarre et al. 2003), 

T Although the quoted Planck sensitivity (AT/T per pixel) is 
slightly better at 30 GHz than at 100 and 143 GHz, the smaller 
pixel size (in arcmin) at high frequencies means that the 100 and 



far exceeding the performances of WMAP, so that an ac- 
curate subtraction of point source contamination at these 
frequencies will be very important. 

A second motivation arises because current models of 
radio-galaxy evolution imply that very young radio sources 
(tens to hundreds of years old) should have radio spectra 
which peak above 5 GHz (the so-called Gigahertz-Peaked 
Spectrum, or GPS sources; O'Dea 1998), with the spectra 
of younger sources peaking at progressively higher frequen- 
cies. Such objects are rare, since they represent a very short- 
lived evolutionary stage, but are predicted to be extremely 
luminous and hence easily detectable out to high redshift in 
high-frequency radio surveys. Until now, models for the lu- 
minosity evolution of young radio galaxies have been almost 
totally unconstrained by the available data; and the need 
for a large, homogeneously-selected sample of GPS sources 
is well-recognized (e.g. Snellen et al. 2000). An analysis of 
evolutionary properties of GPS galaxies has recently been 
carried out by Tinti & De Zotti (2006). 

In addition to these science goals, 95 GHz observations 
are relevant to the calibration strategy for the Atacama 
Large Millimetre Array (ALMA; De Breuck 2005), which 
will need to do rapid switching between object and cali- 
brators to achieve the required phase stability. This ideally 
requires a network of calibrators no more than 1-2 degrees 
apart and stronger than 50-100 mjy at 90 GHz. The ex- 
isting proposal (Holdaway et al. 2004) is to find 90 GHz 
calibrators by selecting candidates from low- frequency (1- 
5 GHz) surveys and using ALMA itself to measure their 
90 GHz flux density. We have already suggested (Sadler et 
al. 2006) that a more efficient alternative would be to select 
ALMA calibrators from the high-frequency Australia Tele- 
scope 20 GHz (AT20G survey). One aim of the present study 
is therefore to measure the surface density of AT20G sources 
at 90-100 GHz. 

Early work in this area dates back to the 1970s. Owen et 
al. (1978) used the NRAO 11m radio telescope at Kitt Peak 
to measure the 90 GHz flux density of 237 strong sources 
which were known to have flat or inverted radio spectra at 
centimetre wavelengths. They also made near-simultaneous 
observations of about 100 of these sources at 1.4, 4.6, 15 and 
22 GHz using dishes of the NRAO VLA. Their main findings 
were that most of their sources had flat spectra (a ~ 0) over 
the whole frequency range 1-90 GHz, and that the level of 
radio variability at 90 GHz was lower than expected (most 
sources varied by less than 30% on a one-year timescale). 
Owen, Helfand & Spangler (1981) observed 25 strong ex- 
tragalactic millimetre sources at ~ 1 keV with the Einstein 
X-ray obsetvatory, and found a strong correlation between 
their X-ray and 90 GHz flux densities. They noted that 
this is consistent with a synchrotron self-Compton model in 
which the X-ray emission arises from the Compton scatter- 
ing of millimetre radio photons off the relativistic electrons 
which created them. 

Steppe et al. (1988) observed a sample of 294 strong 
(> 1 Jy at 5 GHz) extragalactic sources at 90 GHz with the 
IRAM 30 m telescope, and also measured 20-90 GHz spec- 
tral indices for a subset of these objects. They found that 



143 GHz channels are far more sensitive when scaled to the same 
angular resolution. 
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Date 


ATCA 


^ccn 


Notes 


2005 Jul 1-3 


H75A 


18.8, 21.1, 93.5, 95.6 


2x24 hr 


2006 Sep 16-19 


H75A 


18.8, 21.1, 93.5, 95.6 


3xl2hr 



Table 1. Log of our ATCA observations at 95 GHz. 



Run N 20 GHz flux Spectral Notes 

density index cut 

2005 59 > 50mJy af > 'Inverted-spectrum' 

2005 9 > 300mJy ... AT calibrator sample 

2006 70 >150mJy ... Flux-limited sample 



only about 5% of compact, fiat-spectrum radio sources had 
inverted (a > 0.1) spectra at 20-90 GHz, and that most of 
these were variable sources undergoing an outburst. 

More recent continuum work at 90-100 GHz has mainly 
focussed on variability studies of targeted samples of ob- 
jects such as QSOs and blazars (e.g. Landau et al. 1980; 
Terasranta et al. 2004; Hovatta et al. 2007). All these stud- 
ies have been based on sources pre-selected at frequencies 
of 5 GHz or below. The WMAP catalogue (Hinshaw et al. 
2007) provides the first complete sample of radio sources se- 
lected in the 90 GHz band, but only 121 of the 323 discrete 
sources listed in the 3-year catalogue are detected at 94 GHz 
(limiting flux density 0.8 Jy). 

Throughout this paper, we use Ho = 70 km s _1 Mpc -1 , 
Q. m = 0.3 and Q. A = 0.7. 



2 ATCA 95 GHZ OBSERVATIONS 

The data analysed in this paper were obtained at the Aus- 
tralia Telescope Compact Array (ATCA) in two observing 
runs during 2005 and 2006 (see Table 1 for details of the 
observations) . 

The ATCA was used in mosaic mode to minimize the 
slewing time between sources, with the target sources ob- 
served in groups of 5-7 objects which were close (within 
about 10 degrees on the sky) in both right ascension and 
declination. 

2.1 Sample selection 

All our target sources were selected from the Australia Tele- 
scope 20 GHz (AT20G) survey of the southern sky (Ricci 
ct al. 2004; Sadler et al. 2006) as noted in §1. The target 
selection criteria were different in 2005 and 2006, and are 
summarized in Table 2. 

An important feature of the AT20G survey is that 
sources selected at 20 GHz are also observed near- 
simultaneously at 5 and 8 GHz, so that radio spectral in- 
dices measured in the 5-20 GHz frequency range are not 
affected by variability. A recent description of the AT20G 
data pipeline is given by Massardi et al. (2008). 

In July 2005 we observed a sample of 59 sources se- 
lected from the AT20G survey to have inverted (rising) 
radio spectra between 8 and 20 GHz. Nine strong AT20G 
sources which are ATCA calibrators were also observed to 
provide a consistency check on our final flux-density scale 
at 95 GHz. All these sources were also chosen to lie at decli- 
nations between —30° and —50° to ensure that they could 
be observed at low airmass. Our aim was to see how many 
of the 'inverted-spectrum' AT20G sources showed a high- 
frequency turnover above 20 GHz and so were candidates for 
extreme GPS sources. 



Table 2. Sample selection for the 2005 and 2006 runs. All targets 
were also chosen to lie in the declination range —30° to —50° to 
ensure that they could be observed at low airmass. 



In September 2006, we observed a flux-limited (S20 > 
150 mjy) sample of 70 AT20G sources, again between dec- 
lination —30° and —50°. Our aim here was to use the dis- 
tribution of 20-95 GHz radio spectral indices to estimate 
the overall surface density of extragalactic radio sources at 
95 GHz. Eight sources belonged to both the 2005 and 2006 
samples, giving a total of 130 individual sources with 95 GHz 
measurements in the two-year program. 

2.2 Antenna configuration and observing 
frequencies 

We used the ATCA in its most compact configuration, the 
hybrid H75A array, which has a maximum baseline of 75 m 
and is designed to provide good uv coverage and phase sta- 
bility for short observations at millimetre frequencies. The 
ATCA field of view is roughly 2.5arcmin at 20 GHz and 
30 arcsec at 95 GHz, and the angular resolution in the H75 
array is around 40 arcsec at 20 GHz and 10 arcsec at 95 GHz. 

Since one of our goals was to measure reliable 20- 
95 GHz spectral indices for the AT20G sources, and some 
AT20G sources are known to be variable at 20 GHz (Sadler 
et al. 2006), we observed all our target sources at both 20 
and 95 GHz in the 2005 and 2006 runs. The observing fre- 
quencies (IFs) listed in Table 1 were selected to lie in regions 
where the ATCA receivers have their best Tsys performance. 

2.3 Calibration techniques 

At the time of these observations the ATCA 3 mm system 
did not have a noise diode, so to place visibility amplitudes 
on a temperature scale a vane (paddle) needed to be inter- 
posed in front of the receiver in the antenna vertex room at 
regular intervals for Tsys calibration. This was done at the 
start of each observing session, and also repeated just before 
starting the observation of each source group (at intervals 
of about half an hour). 

2.3.1 Flux-density calibration 

At 20 GHz, the absolute flux-density scale was set using 
the standard ATCA primary calibrator PKS 1934-638. At 
95 GHz, however, PKS 1934—638 is far too weak to provide 
a reliable calibration of the flux-density scale. We therefore 
followed the usual practice and used the planet Uranus as 
our flux calibrator at 95 GHz. Although the angular size of 
this planet is significantly smaller than the ATCA beam, 
it cannot be assumed to be a point source. A model of its 
intensity as a function of frequency and baseline length is 
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Figure 2. The effect of antenna pointing errors on the measured 
flux density at 20 and 95 GHz. For our 2005 and 2006 runs, point- 
ing and tracking errors arc expected to be less than 4.0 arcsec 
(shown by the dotted vertical line), except on rare occasions when 
the wind is both strong and gusty. Pointing and tracking errors 
should therefore have at most a small effect on the measured flux 
densities and 20-95 GHz spectral index distribution. 



therefore used to perform the calibration. To minimize any 
gain-elevation effects (i.e. possible variations in antenna effi- 
ciency as a function of elevation) on the primary calibration, 
we always observed the planet close to transit and followed 
or preceded by a target source group at a similar elevation. 



2.3.2 Antenna pointing calibration 

The ATCA primary beam at 95 GHz is only 30 arcsec 
FWHM, and the pointing and tracking requirements for op- 
eration at this high frequency have been discussed in de- 
tail by Kesteven (1998). The effects of a pointing offset can 
be seen in Figure 2. The main factors affecting the posi- 
tional accuracy of the ATCA are pointing and tracking er- 
rors (typically 1.5 arcsec rms on each axis when reference 
pointing is used), atmospheric fluctuations (up to 3.0 arcsec 
rms for sources observed close to the zenith) , and wind load- 
ing (which starts to have a significant effect when the wind 
speed rises above 10-15 kmhr - , especially when conditions 
are gusty). Also contributing to the error budget is the un- 
certainty in the individual AT20G source positions (typically 
less than 1.0 arcsec). 

The uncorrected ATCA antenna pointing errors in typ- 
ical observing conditions can be up to 20 arcsec. It is there- 
fore essential to use reference pointing when observing at 
high frequencies. A reference pointing check using a strong 
phase calibrator was therefore done whenever moving to a 
new source group. This allowed a new pointing solution to 
be calculated and applied to the sources in each group, re- 
ducing the final pointing errors to a few arcsec or less. We 
generally did our reference pointing calibrations at 20 GHz 
rather than 95 GHz, because the higher signal-to-noise ratio 
and larger field of view at 20 GHz allows a more accurate 
solution. 



2.4 2005 observing run 

Both the effective detection limit for continuum sources at 
95-100 GHz and the accuracy of the derived flux densities 
are strongly affected by atmospheric transparency, phase 
stability and decorrelation, and the general performance of 
the ATCA 3 mm system for weak continuum sources was 
unknown at the time of our 2005 run. We therefore took a 
conservative approach and repeated the same set of obser- 
vations on each of the two days we were allocated. 

Since all our target sources were known (from AT20G 
images at 20 GHz) to be point sources at high frequency, 
we chose to measure flux densities using a triple correlation 
technique rather than imaging our target fields. The triple 
correlation technique is discussed briefly in §3 of this paper, 
and in more detail in a companion paper by Ricci, Sault & 
Ekers (in preparation). 

In the 2005 run, each group of target sources was ob- 
served three times at different hour angles to provide a range 
in uv coverage. The total integration time for each source 
at 95 GHz was 3x4 minutes on each day. On the second 
day, each group of sources was observed at both 95 GHz 
and 20 GHz, with the total observing time at 20 GHz being 
3x1 minute. To minimize the effects of airmass, all obser- 
vations were made at hour angles less than ±1.5 hours from 
the zenith. 

2.5 2006 observing run 

During our 2005 run, we found that observing conditions at 
95 GHz were significantly better at night than during the 
day. This appears to be generally the case at the ATCA, 
even during the winter months, as discussed recently by Mid- 
delberg et al. (2006). We therefore chose to carry out our 
2006 observing campaign in three consecutive 12 hr night- 
time slots. The array configuration and frequency set-up 
were the same as in 2005. 

Our 2006 sample of 70 AT20G sources was split into 
nine groups which were observed at 20 and 95 GHz on either 
the first or the second nights. Each source was observed be- 
tween two and four times, with the integration time for each 
observation being 80 s at 95 GHz and 40 s at 20 GHz. The in- 
tegration times were shorter than in 2005 because most of 
the 2006 sources were significantly brighter at 20 GHz than 
those observed in 2005. On the third night, we reobserved 32 
sources which did not appear to show an obvious 95 GHz de- 
tection on the first two days. These objects were re-observed 
for 3min cuts at 95 GHz. 



3 DATA REDUCTION 

The raw visibility data were reduced using the astronomical 
software package MIRIAD (Sault et al. 1995). The 95 GHz 
and 20 GHz bandpasses were calibrated using the strong 
ATCA high-frequency calibrator PKS 1921-293. 

3.1 Flux— density measurements 

Atmospheric turbulence creates strong phase instability in 
the 95 GHz band and target sources and phase calibrators 
need to be switched very frequently (i.e. on timescales of a 



Radio sources at 95 GHz 5 



Flux— density error estimates 
20 GHz 95 GHz 



2005 5.4% 11.8% 

2006 10.7% 10.2% 



Table 3. Error estimates for the flux densities measured in 2005 
and 2006, derived using the methods of Ricci et al. (2007). 



few minutes or less) to avoid the loss of total intensity flux 
through phase decorrclation. 

Because of the large number of sources we wanted to 
observe in both 2005 and 2006, it would have been difficult 
to spend enough time on phase calibrators to correct for 
phase decorrelation in the standard way. We therefore chose 
not to observe any secondary (phase) calibrators in our 2006 
run. Instead, we chose to reduce our 2005 and 2006 data 
using the triple correlation (phase closure) method as noted 
earlier. 

The triple correlation method has also been used by 
Sault et al. when monitoring flux densities of ATCA calibra- 
tors at 3 mm, and is robust against phase variations during 
the target observation. It is described in detail by Cornwell 
(1987) and Thompson et al. (2001). The application of this 
technique to ATCA millimetre data is discussed by Ricci et 
al. (in preparation), who have also used a set of Monte Carlo 
simulations to test the robustness of this technique for the 
detectability and flux density measurements of the 95 GHz 
data presented in this paper. 

Table 3 lists the estimated flux-density errors for our 
2005 and 2006 data. These were calculated from the scatter 
in the triple-correlation amplitudes by taking different scans 
and different spectral windows (IFs) of the same source sep- 
arately. At 20 GHz, the flux density errors are significantly 
higher in 2006 than in 2005. This is probably because the 
2005 20 GHz data were taken in very good observing con- 
ditions, whereas the 2006 data at the same frequency were 
much noisier. 



3.2 Data tables 

Tables 4 and 5 list the individual sources observed at 20 and 
95 GHz in our 2005 and 2006 runs. The main information 
contained in these tables is as follows: 

(i) The AT20G working name, used to track individual 
sources through the data pipeline. A letter 'C after the name 
indicates that the source is also an ATCA calibrator. 

(ii) The J2000 radio position of the source, measured from 
AT20G 20 GHz images. 

(iii) Alternate source name, where available. 

(iv) Blue (Bj) magnitude of the optical counterpart if 
any, from the Supercosmos catalogue (Hambly et al. 2001). 
As discussed in §4.7, an optical object is accepted as the 
radio ID if it lies within 2.5 arcsec of the AT20G radio posi- 
tion. 

(v) Supercosmos classification of the optical ID, where 1 
indicates a galaxy and 2 a stellar object (candidate QSO). 

(vi) Simultaneous ATCA 95 GHz and 20 GHz flux densi- 
ties measured in 2005/6 in this project. 



(vii) Near-simultaneous ATCA 20, 8 and 5 GHz flux den- 
sities measured by the AT20G team in November 2004. The 
listed values are from the July 2007 run of the AT20G data 
pipeline. 

(viii) Lower-frequency flux-density measurements from 
the 1.4 GHz NVSS (Condon et al. 1998) and 843 MHz 
SUMSS (Mauch et al. 2003) catalogues. The SUMSS val- 
ues are taken from v2.0 of the catalogue (released August 
2007). 

(ix) The spectroscopic redshift of the optical counterpart, 
where this is available either from the NASA Extragalactic 
Database (NED) or from the 6dFGS Third Data Release 
(Jones et al., in preparation). 

(x) NED classification as QSO (Q) or galaxy (G), where 
available. These classifications are taken from papers in the 
literature and usually (but not always) agree with the Su- 
percosmos classification. 

Appendix A contains some notes on the properties of 
individual sources. 



4 DATA ANALYSIS 

4.1 Stability of the 95 GHz flux-density 
measurements 

A comparison of the 95 GHz flux densities measured on the 
first and second days of our 2005 run is shown in Figure 3. 
The flux densities measured on the two days agree to within 
10% (i.e. within the quoted error bars from Table 3), even 
though conditions on the first day were cloudy with poor 
phase stability. We therefore conclude that our techniques 
for measuring 95 GHz flux densities yield results which are 
stable and reproducible. This is important, because we are 
measuring continuum sources which are significantly weaker 
than have previously been observed with the ATCA in the 
95 GHz band. 



4.2 The inverted— spectrum and flux— limited 
samples 

As noted in Table 2, our 2005 observations included 59 
inverted-spectrum sources which were selected to have ris- 
ing radio spectral index (a > 0) between 8 and 20 GHz. The 
sample selection was done in July 2005, using preliminary 
flux-density measurements from a reduction of the AT20G 
data carried out in late 2004. 

After our 2005 observations were completed, the AT20G 
team implemented a new data pipeline which included an 
improved calibration process and a stringent set of checks 
designed to identify poor-quality data. As a result of these 
checks, 13 galaxies in our 2005 sample (all in the RA range 
17-22 h) had their 8 or 20 GHz AT20G flux measurements 
flagged as unreliable. Since we could no longer be sure that 
these sources qualified as 'inverted spectrum', they were not 
included in our analysis. Two other sources no longer had 
rising 8-20 GHz spectra after recalibration, and were also 
removed. 

Finally, we tightened the selection criteria for the in- 
verted spectrum sample to require that the sources had 
a > at 5-8 GHz as well as 8-20 GHz. The aim was to 
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Name — S20 (mjy) — S95 (mjy) - — Spectral index a|° - 







2004 


± 


2005 


± 


2006 


± 


2005 


± 


2006 


± 


2005 


± 


2006 


± 


0013- 


-3954 


1609 


33 


1420 


77 


984 


105 


716 


84 


496 


51 


-0.44 


0.12 


-0.44 


0.14 


0106- 


-4034 


2146 


79 


2632 


142 


4160 


445 


1775 


209 


2199 


224 


-0.25 


0.12 


-0.41 


0.14 


0136- 


-4044 


350 


11 


380 


21 


274 


29 


344 


41 


190 


19 


-0.06 


0.12 


-0.23 


0.14 


0421- 


-4817 


452 


13 


513 


28 


321 


34 


300 


35 


212 


22 


-0.34 


0.12 


-0.27 


0.14 


2056- 


-4714 


1711 


55 


1993 


108 


3246 


347 


1086 


128 


2453 


250 


-0.39 


0.12 


-0.18 


0.14 


2151- 


-3027 


1846 


58 


2108 


114 


1402 


150 


688 


81 


849 


87 


-0.72 


0.12 


-0.32 


0.14 


2158- 


-3013 


350 


10 


388 


21 


844 


90 


327 


39 


903 


92 


-0.11 


0.12 


+0.04 


0.14 


2235- 


-4835 


1994 


75 






1775 


190 


2605 


307 


1076 


110 






-0.32 


0.14 



Table 6. Flux-density and spectral-index measurements for the eight sources from Tables 4 and 5 which were observed in both 
the 2005 and 2006 runs. 2235-4835 was not observed at 20 GHz in 2005. 
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Figure 3. Comparison of two independent 95 GHz measurements 
(taken one day apart) of AT20G sources observed in our July 2005 
run. Dotted lines show a 10% deviation from equality. 



restrict the sample to objects whose radio spectra were gen- 
uinely inverted over the full 5-20 GHz AT20G frequency 
range, and to avoid contamination by sources whose spectra 
were fiat (a ~ 0) at 5-20 GHz but which happened to have a 
slightly higher flux measurement at 20 GHz than at 8 GHz. 
A further seven (mostly weak) sources were removed by this 
test. All 59 sources observed in 2005 are listed in Table 4, 
with the 39 sources included in our final inverted-spectrum 
analysis clearly identified. 

Table 5 lists the 70 sources in our 2006 'flux-limited' 
sample, which were selected to have 20 GHz flux densities 
above 150 mjy. One object (2257—3627, the nearby galaxy 
IC 1459) later had its 2004 AT20G flux density at 20 GHz 
flagged as unreliable. Although it is listed in Table 5, it was 
not included in our analysis. Our final flux-limited sample 
therefore contains 69 objects. 



Data set 


Freq 


Time interval 


N 


Median var. 




GHz 






index (%) 


AT20G Pilot 


20 


2003 - 2004 


108 


6.9 


Inverted-spec 


20 


2004 - 2005 


37 


6.4 


Flux-lim. 


20 


2004 - 2006 


69 


6.9 


05-06 repeats 


20 


2004 - 2006 


8 


20.1 


05-06 repeats 


20 


2005 - 2006 


7 


20.3 


05-06 repeats 


95 


2005 - 2006 


8 


20.2 



Table 7. Median variability index at 20 and 95 GHz for several 
sets of data observed with the ATCA. The AT20G Pilot results 
for 2003-4 are taken from Sadler ct al. (2006). The small sample 
of 'repeat' sources which we observed in both 2005 and 2006 have 
a significantly higher median variability index than either of the 
parent samples from which they came. 



4.3 Variability at 20 and 95 GHz 

Sadler et al. (2006) examined the variability of a flux-limited 
(S20 > 100 mjy) sample of 108 AT20G sources by comparing 
measurements from 2003 and 2004. They found that the 
general level of variability at 20 GHz was quite low, with a 
median variability index for their sources of 6.9%. Only five 
sources varied by more than 30% over the one-year interval. 

We can also use the 2004 AT20G measurements as a 
benchmark to measure the variability of the sources in Ta- 
bles 4 and 5 at 20 GHz. To do this, we calculated a debiased 
variability index for each object as described in §6 of Sadler 
et al. (2006). The results are summarized in Table 7. For 
both the 'inverted-spectrum' sample observed in 2005 (Ta- 
ble 4) and the flux-limited sample observed in 2006 (Table 
5), the median variability index at 20 GHz is close to that 
found by Sadler et al. (2006) for their flux-limited source 
sample. 

In contrast to the results at 20 GHz, we have very little 
data so far on variability at 95 GHz. There are only eight 
sources in common between our 2005 and 2006 data sets 
(see Table 6), and these turn out by chance to be drawn 
mainly from the more variable objects in their parent sam- 
ples. Any results based on this small sample are no more 
than suggestive, but we see no evidence at this stage that 
our sources are more variable at 95 GHz than at 20 GHz (see 
Table 7). The spectral index values in Table 6 also suggest 
that the measured 20-95 GHz spectral index may not change 
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ATCA, 95 GHz flux density (Jy) 




-1.5 -1 -0.5 

20—95 GHz spectral index a 



Figure 4. Comparison of our measured ATCA 95 GHz flux den- 
sities with those listed in the WMAP 3-year source catalogue 
(Hinshaw ct al. 2007) for sources in common (see also Table 8). 
Filled and open circles show data from the 2005 and 2006 ATCA 
runs respectively. 



Figure 5. Histogram of 20— 95 GHz spectral index a for (a) the 
inverted— spectrum sample observed in 2005 and (b) the flux- 
limited sample observed in 2006. 



significantly when a sources varies in flux density at these 
frequencies, though the size of the error bars makes it dif- 
ficult to draw any strong conclusion from a sample of only 
seven repeat observations. 



4.4 Comparison with the WMAP 3— year Source 
Catalogue 

Hinshaw et al. (2007) have recently released a catalogue of 
323 radio sources detected by the WMAP satellite (Ben- 
nett et al. 2003) in five frequency bands at 23, 33, 41, 61 
and 94 GHz. The WMAP catalogue identifies bright sources 
independent of their presence in other surveys at lower fre- 
quency. When a source is identified with > 5a confidence 
in any of the five bands, the flux densities in other bands 
are also listed if they are > 2a. Only 121 of the 323 cata- 
logued WMAP sources have flux densities listed at 94 GHz, 
and almost all of these are 23 GHz detections. 

Table 8 lists the thirteen sources in Tables 4 and 5 which 
were also detected by WMAP at 94 GHz. These objects 
are also plotted in Figure 4. The mean flux-density ratio 
(Satca/Swmap) for the objects in Table 8 is 1.18 ± 0.14 
at 20/23 GHz and 1.16 ± 0.06 at 95/94 GHz. Although the 
comparison sample is small and the sources were measured 
at different epochs, the results strongly suggest that the 
ATCA and WMAP data are on a common flux-density scale 
to within about 15% at 95 GHz. This also appears to con- 
firm that our 95 GHz ATCA flux measurements have not 
been strongly affected by pointing and tracking effects (see 
§2.3.2). 



4.5 The 20—95 GHz spectral— index distribution 

As a first step in investigating the high-frequency properties 
of our sample, we calculated the two-point radio spectral 
index a|o for each of the objects in Tables 4 and 5. For the 
small number of 20 GHz sources which were undetected at 
95 GHz, we used the upper limits in flux density to calculate 
a limiting value of a. Figure 5 shows the distribution of a§o 
for the two samples. Only two sources in the flux-limited 
sample (0227—3026 and 2158—3013) show a rising spectrum 
from 20-95 GHz. Even in the inverted-spectrum sample (all 
of which were selected to have a|° > 0), only seven of the 
37 sources have a^o > 0. This contrasts strongly with the 
spectral-index distribution of AT20G sources at 8-20 GHz, 
where 40% have rising spectra with q|° > (Sadler et al. 
2006). As discussed by Lowe et al. (2007) this may be a 
selection effect, since the number of rising-spectrum sources 
in a sample selected at frequency v will usually be greater 
when the spectral index is measured between v and a lower 
frequency v\ than when it is measured between v and a 
higher frequency v%. 

The median 20-95 GHz spectral index for sources in our 
flux-limited sample is —0.39 and, as can be seen from Fig- 
ure 6, this does not change significantly with 20 GHz flux 
density over the range probed by our ATCA observations 
(S20 > 150 mjy). As we will show in §5, this allows us to 
use the 20 GHz extragalactic source counts measured from 
the AT20G survey, together with the observed distribution 
of alo) to derive the overall radio source counts at 95 GHz. 

We note that our median 20-95 GHz spectral index of 
—0.39 is much flatter than the value of —0.89 measured at 
15-43 GHz by Waldram et al. (2007). The Waldram et al. 
(2007) sources are much fainter than those in our sample, 
since most of them have 22 GHz flux densities below the 
AT20G survey limit of 40 mjy. The work of Waldram et al. 
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— 2006 - - 2005 - — 2004 — WMAP 



AT20G name 


Alt. name 


S95 


± 


S20 


± 


S95 


± 


S20 


± 


S20 


± 


S94 


± 


S23 


± 






Jy 




Jy 




Jy 




Jy 




Jy 




Jy 




Jy 




0012-3954 


WMAP 202 


0.50 


0.05 


0.98 


0.10 


0.72 


0.08 


1.42 


0.08 


1.61 


0.03 






1.4 


0.06 


0026—3512 


PMNJ0026— 3512 










1.12 


0.13 


1.27 


0.07 


1.12 


0.02 






1.0 


0.08 


0106-4034 


WMAP 171 


2.20 


0.22 


4.16 


0.45 


1.78 


0.21 


2.63 


0.14 


2.15 


0.08 


1.3 


0.3 


1.8 


0.05 


0334-4008 


WMAP 146 


0.35 


0.04 


0.55 


0.06 










1.27 


0.04 






1.6 


0.06 


0440-4333 


WMAP 147 










0.44 


0.04 


1.53 


0.16 


1.95 


0.05 






3.0 


0.07 


0455-4615 


WMAP 151 










3.07 


0.36 


4.42 


0.24 


4.16 


0.12 


1.8 


0.4 


3.7 


0.06 


0538-4405 


WMAP 148 










8.08 


0.95 


7.33 


0.40 


5.29 


0.22 


4.6 


0.4 


5.4 


0.06 


1102-4404 


PKS1059-438 










0.49 


0.06 


0.49 


0.03 


0.77 


0.02 






0.6 


0.05 


1107-4449 


WMAP 166 










0.65 


0.08 






1.67 


0.04 


1.0 


0.4 


1.6 


0.05 


1427-3305 


WMAP 193 










1.47 


0.17 


0.87 


0.05 


0.77 


0.02 






3.2 


0.07 


1937-3958 


PKS1933-40 










1.29 


0.15 


1.67 


0.09 


1.76 


0.05 






0.7 


0.08 


2056-4714 


WMAP 208 


2.45 


0.25 


3.25 


0.35 


1.09 


0.13 


1.99 


0.11 


1.71 


0.06 


2.4 


0.9 


1.7 


0.06 


2235-4835 


WMAP 206 


1.08 


0.11 


1.78 


0.19 


2.61 


0.31 






1.99 


0.08 


2.1 


0.5 


1.7 


0.06 



Table 8. Comparison of flux densities for sources in common with the WMAP 3-year catalogue (Hinshaw ct al. 2007). The WMAP 
data were taken over the period 2001 Aug to 2004 Aug inclusive and so pre-date the ATCA observations, which were made in 2004 
Oct, 2005 Jul and 2006 Sep. 
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Figure 6. Plot of 20-95 GHz spectral index a|jj against 20 GHz 
flux density for the flux— limited sample of sources in Table 4. The 
horizontal dashed line shows the median spectral index of —0.39, 
while the vertical dotted line shows the flux limit of 150 mjy 
used to select the sample. A few sources lie to the left of this line 
because their 20 GHz flux density has decreased since the AT20G 
observations were made in 2004. 



(2007) is therefore complementary to our study, and suggests 
that the 20-95 GHz spectral-index distribution of AT20G 
sources may steepen significantly at flux densities below 50- 
100 mjy. 



4.6 Properties of the 5—95 GHz radio spectra 

For the sources in this study, we have near-simultaneous 
data at 5, 8 and 20 GHz from the AT20G survey as well as 
simultaneous data (at a different epoch) at 20 and 95 GHz 
from our own ATCA observations. The mean and median 
values of the spectral index in each frequency interval are 
listed in Table 9. Both samples show a general steepening 
of the radio spectrum with increasing frequency. A similar 



-2 



Upturn above 20 GHz 



Inverted (rising) 




Steep (falling) 



Spectral peak (GPS) 



-2 -1 1 2 

Spectral index (B.6 to 20 GHz) 

Figure 7. Two-colour radio spectral-index plot for the flux- 
limited (filled circles) and inverted-spectrum (open circles) 
AT20G samples. Only three sources still have a rising spectrum 
from 20-95 GHz, and almost all the sources with inverted spectra 
at 8-20 GHz have peaked and turned over by 95 GHz. 



effect has been noted by Waldram et al. (2007) and Massardi 
et al. (2008). 

Figure 7 shows a radio two-colour diagram comparing 
the 8-20 and 20-95 GHz spectral indices. As discussed by 
Sadler et al. (2006), this diagram provides a useful way of 
classifying objects which (like many of the AT20G sources) 
have curved, rather than power-law radio spectra. The di- 
agonal dotted line in Figure 7 shows the power-law relation. 
As noted earlier, inverted-spectrum objects (which account 
for 18% of all AT20G sources at 8-20 GHz), are almost ab- 
sent (< 2%) at 20-95 GHz. 
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Sample 


N 


Gal 


QSO 


Faint/blank 


Median Bj 


N z 


Median z 












(mag) 






Inverted spectrum 


34 


7 (21 ± 8%) 


23 (68 ± 14%) 


4 (11 ± 7%) 


19.84 


12 (35%) 


0.56 


Flux-limited 


69 


11 (16 ±5%) 


50 (73 ± 10%) 


8 (11 ±4%) 


19.44 


32 (46%) 


1.11 



Table 10. Optical properties of the two samples (excluding sources for which no optical ID was attempted). N z is the number 
of sources in each sample which have a published optical redshift. Objects classed as 'Faint/blank' have optical IDs fainter than 
Bj=22.0. 



Frequency Median Mean ± n Notes 

(GHz) a a 

5-8 -0.11 -0.17 0.06 62 This paper 

8-20 -0.24 -0.25 0.05 62 This paper 

20-95 -0.39 -0.45 0.04 69 This paper 



Table 9. Mean and median values of the radio spectral-index 
distribution for our 2006 flux-limited sample calculated for dif- 
ferent frequency ranges. The three sources undetected at 95 GHz 
have been removed from the calculation of the mean value at 
20-95 GHz. 



4.7 Optical properties of the 95 GHz sources 

We used the Supercosmos catalogue (Hambly et al. 2001) to 
obtain optical identifications for the radio sources observed 
in this study. We accepted an optical object as the correct ID 
if it was brighter than Bj = 22 mag and lay within 2.5 arcsec 
of the radio position, since Monte Carlo tests imply that at 
least 97% of such objects are likely to be genuine associ- 
ations (Sadler et al. 2006). Optical identifications were not 
attempted for three sources within 10 degrees of the Galactic 
plane, and one source in the nux-limted sample (0320—3837) 
was so close to a bright foreground star that no optical iden- 
tification was possible. 

Table 10 summarizes the optical properties of the two 
samples. While QSOs make up the majority of the optical 
IDs in both samples, the inverted-spectrum sample has a 
slightly higher fraction of galaxies and a median value of Bj 
0.4 mag fainter than the flux-limited sample. 

Published redshifts are available for fewer than half the 
optical IDs, so the redshift distribution needs to be inter- 
preted with caution. However, the lower median redshift 
for the inverted-spectrum sample is consistent with it con- 
taining a higher fraction of galaxies (which generally lie at 
z < 0.5 if visible on DSS images) than the flux-limited sam- 
ple. 

4.8 A correlation between 95 GHz flux density 
and optical magnitude? 

Owen & Mufson (1977) found a correlation between optical 
magnitude and millimetre (90 GHz) flux density for a sample 
of flat-spectrum QSOs, and noted that this was surprising 
because of the 20-25 year interval between their 90 GHz mea- 
surements and the sky survey plates from which the optical 
magnitudes were measured. They showed that the correla- 
tion was significant at the 99% level at 90 GHz, yet the 5 GHz 
flux densities of the same sources were uncorrelated with op- 
tical magnitude. They interpreted this as evidence that the 
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Figure 8. Plot of optical Bj magnitude from the SuperCOS- 
MOS catalogue against 95 GHz flux density for objects in our 
flux-limited sample. Sources identified with QSOs are shown by 
stars, and galaxies by open circles. Upper limits are indicated by 
open triangles, and the horizontal line at 50 mjy shows the typical 
detection limit for our 2006 observations. 



sources were optically thin at millimetre wavelengths but 
optically thick at centimetre wavelengths, and that a cor- 
relation existed between the radiation mechanisms in the 
optical and millimetre regimes. 

Figure 8 shows a plot of 95 GHz flux density against 
optical Bj magnitude for our flux-limited sample, with sep- 
arate symbols used for QSOs and galaxies. Blank fields are 
indicated by upper limits in Bj. We see a trend similar to 
that found by Owen & Mufson (1977), in the sense that 
brighter 95 GHz sources are associated with brighter optical 
objects. The correlation for QSOs in Figure 8 is statistically 
significant at a level of > 99.6% in a rank correlation test. 

One possible explanation is that many of the strongest 
95 GHz sources are brightened by relativistic beaming at 
both optical and millimetre wavelengths, which would re- 
quire that they are viewed at very small angles to the jet 
axis. In this orientation, the millimetre continuum emis- 
sion can be Doppler-boosted by factors of up to twenty 
(Lahtecnmaki & Valtaoja 1999) and the optical continuum 
brightened by at least one magnitude (Browne & Wright 
1985). 
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In principle, optical spectroscopy can provide a test 
of the relativistic beaming model, since QSOs which are 
strongly beamed should have a lower line to continuum ratio 
than less beamed objects (Baker & Hunstead 1995). This 
would require a more complete and uniform set of optical 
spectra than is currently available for our sample. 

4.9 Sources with extended low— frequency radio 
emission 

SUMSS images at 843 MHz are available for all our 95 GHz 
sources, and 1.4 GHz NVSS images for those north of dec- 
lination -40°. The SUMSS and NVSS images have similar 
resolution (45 arcsec) and sensitivity (detection limits of a 
few mjy) . Although all the sources in Tables 4 and 5 are un- 
resolved at 20 and 95 GHz, several of them show extended 
emission on scales of 45-100 arcsec in the low-frequency 
SUMSS and NVSS images (see Appendix A). From the 
NVSS/SUMSS images alone it is unclear whether the ex- 
tended low-frequency emission arises only from the AT20G 
object, or includes a contribution from a nearby confusing 
source which is below the detection threshold at high fre- 
quencies. 

Interestingly all the sources with extended low- 
frequency emission are members of the inverted-spectrum 
sample. Five of the 37 sources in the inverted-spectrum sam- 
ple (14%) have extended emission in the SUMSS images, 
whereas none of the 70 sources in the flux-limited sample 
do. The reason for this difference is unclear at present. It 
may simply reflect the lower median redshift of the inverted- 
spectrum sample (see Table 10), but could also hint at 
genuine physical differences between the inverted-spectrum 
sources (and/or their environment) and the overall high- 
frequency radio population. 



5 ATCA 95 GHZ SOURCE COUNTS AND 
COMPARISON WITH EARLIER 
PREDICTIONS 

As noted earlier, there are several reasons for wanting to 
measure the surface density of extragalactic radio sources at 
95 GHz apart from the study of the intrinsic physical prop- 
erties of these objects. In an ideal case one would perform 
an all-sky survey at 90-100 GHz, but it is currently difficult 
to reach faint flux levels with such a survey. As noted in §1, 
the primary beam of most radio interferometers and single- 
dish telescopes operating at mm-wavelength is tiny, making 
the mosaicing of any large area of sky an incredibly time- 
consuming task't . Additional challenges include the high sys- 
tem temperatures at 90-100 GHz and changing weather con- 

t For example, a blind survey carried at 100 GHz with the same 
observational setup used for the AT20G survey at 20 GHz (3 base- 
lines; 8 GHz analogue correlator and fast meridian scanning; see 
e.g. Sadler ct al. 2006) on the same 1500 deg 2 sky strip explored 
by the AT20G Pilot Survey (-70° < 8 < -60°; see Ricci et al. 
2004) would require a month of observing time at the ATCA in- 
stead of the five days needed at 20 GHz. With a detection limit of 
~ 300 mjy, this 100 GHz survey would detect less than one source 
a day. 



ditions which can make the sensitivity level of ground-based 
surveys uneven. 

The Planck satellite scheduled for launch in 2008 will 
survey the whole sky with an estimated 5a point-source 
detection limit of ~ 500 mjy at 30 GHz and ~ 280 mjy at 
100 GHz (Lopez-Caniego et al. 2006). These limits are set by 
confusion, and the rms noise per pixel is much lower (14mJy 
at 100 GHz and 10.2 mjy at 143 GHz (Lamarre et al. 2003). 
As a result, there is a lot of astrophysical information in 
Planck maps below the confusion limit, and some of this can 
be extracted, e.g. using stacking techniques, by making use 
of the AT20G survey and follow-up observations at higher 
frequencies. 

Since a blind survey of the whole sky to low flux levels 
at frequencies above 20 GHz is currently impractical, we use 
a method first proposed by Kellermann (1964) to predict the 
extragalactic source counts at 95 GHz, based on our knowl- 
edge of the observed 20 GHz source counts from the AT20G 
survey and the observed 20-95 GHz spectral-index distribu- 
tion presented in §4.5 of this paper. 



5.1 Methodology 

The technique developed by Kellermann (1964) is based on 
the idea that if one knows both the source counts at one fre- 
quency and the distribution of flux densities for a complete 
sub-sample at a second frequency, it is possible to predict 
the source counts at the second frequency. 

The method makes no assumption about the actual 
spectral shape between the two frequencies, but it is con- 
venient to characterize the flux ratio by the spectral index 
for an assumed power law. It does assume that the distri- 
bution of flux ratios (i.e. the spectral-index distribution) 
correctly characterizes all the source populations present. If 
the spectral-index distribution is a function of flux density, 
or if some class of sources becomes important at the second 
frequency but is not represented at the first frequency, the 
prediction will fail. For this reason it is important to start 
with reliable source counts at a frequency as close as possible 
to the second frequency. 

The method we followed (outlined by Kellermann et al. 
1968) first assumes that the ratio of the flux densities Si 
and S2 at two frequencies v\ and V2 is represented by the 
spectral index a defined as 

s 2 = s lX c^r (i) 

V\ 

The observed spectral index distribution P(a) and the 
counts n(S\,v\) at any given frequency v\ above the flux 
density limit Si can then be used to predict the source 
counts n(S2, vi) at another frequency V2 above the flux den- 
sity limit S2 using the formula: 

/+00 
n[(^)- Q S a> i/i] P(a) da (2) 
"1 

provided that the P{a) is independent of flux density. 

For our flux-limited sample with S20 > 150 mjy and a 
median 20-95 GHz spectral index of —0.39, the equivalent 
flux-density limit at 95 GHz is 82 mjy. 
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Parameter Fitted value Units 



a. 


1 Q2+003 




p 


r, 7a +0.09 
2 - 78 -0.07 






1 09 +0 13 


Jy 


N, 


28.9te 3 8 7 


Jy sr _1 



Table 11. Best— fit parameters for the 20 GHz source counts mea- 
sured from the AT20G survey at declination -90° < 8 < -15° 
(Ricci ct al., in preparation). 



Parameter 


Fitted value 


Units 


a 


1.92 






2.77 




s* 


0.80 


Jy 


N» 


19.45 


Jy sr _1 



Table 12. Parameterization of our predicted 95 GHz source 
counts plotted in Figure 9. The fit was done in the same way 
as for the 20 GHz counts in Table 11. 



5.2 Completeness of the 95 GHz source 
populations 

As noted above, the method we are using to predict the 
95 GHz source counts will give unreliable results if any class 
of sources which is present at 95 GHz is not represented in 
our 20 GHz source counts. 

In our case, any new population of sources appearing 
at 95 GHz would have to have a very rapidly-rising spec- 
trum between 20 and 95 GHz in order to be weaker than 
the AT20G detection limit of 40mJy at 20 GHz. At the 
1 Jy level at 95 GHz, this would require a population with a 
high-frequency radio spectral index a > +2.0. Even at the 
100 mjy level, a high-frequency spectral index a > +0.6 
would be required. 

No sources with such extreme 20-95 GHz spectral in- 
dices are seen in our current 95 GHz data set, apart from a 
small population of optically-thick thermal sources in our 
own Galaxy (which all lie at low Galactic latitude). As can 
be seen from Figures 5 and 7, almost all AT20G sources with 
a > at 8-20 GHz have lower flux densities at 95 GHz than 
at 20 GHz. Bolton et al. (2004) also found no sources with 
a > +0.6 at 15-43 GHz in their sample of sources selected 
at 15 GHz, with flux densities as low as 25 mjy. 

At this stage, therefore, we see no evidence that that our 
20 GHz sample is missing any source population which con- 
tributes to the 95 GHz source counts at levels above about 
100 mjy. 

5.3 Differential radio— source counts at 95 GHz 

Since the 20 GHz counts show some curvature, we made our 
final estimate by fitting the AT20G source counts with a 
double power law of the form: 

N(S) = 2N»/[(S/S,r + (S/S»)' 3 ] (3) 

where the best-fit parameters for the two slopes a and (3, 
the flux density of the break S* and number count at the 
break N, are given in Tables 11 and 12. 

We then derived the differential source counts at 95 GHz 
using the methods outlined in §5.1. The results are shown 
in Figures 9 and 10, and are valid for flux densities above 
80 mjy at 95 GHz. 

5.4 Comparison with previously— published values 

Figures 9 and 10 also compare our 95 GHz differential source 
counts with previous predictions by Holdaway et al. (1994), 
De Zotti et al. (2005) and Waldram et al. (2007). 



Holdaway et al. (1994) observed a sample of 367 flat- 
spectrum radio sources at 90 GHz with the NRAO 12 m tele- 
scope in 1993. Just under 80% of their sample was detected 
above a 3a limit of 90 mjy. Holdaway et al. (1994) then used 
VLA 8.4 GHz measurements of the same objects (Patnaik et 
al. 1992) to calculate the 8-90 GHz spectral index distribu- 
tion for their sample, and combined this with the 5 GHz 
source counts measured by Condon (1984) to estimate the 
source counts at 90 GHz. 

The 8 and 90 GHz observations were separated in time 
by almost four years (the data were taken in 1990 Feb and 
1993 Nov respectively) , and so the measured 8-90 GHz spec- 
tral index distribution is likely to be broadened by source 
variability over the intervening period of almost four years. 
The Holdaway et al. (1994) analysis also depends critically 
on the assumed ratio of core to total flux density at 5 GHz, 
since the 5 GHz source counts measure total flux density but 
their 8-90 GHz data essentially measures the spectral index 
of a compact core. 

Their predicted source counts from Holdaway et al. 
(1994) are valid for the flux-density range S90 > 0.1 Jy, set 
by the detection limit of the high-frequency survey they 
used. 

Radio source counts at 94 GHz have been estimated us- 
ing the evolutionary models described by De Zotti et al. 
(2005), which take into account the epoch-dependent ra- 
dio luminosity functions of the various source populations 
which contribute to the high-frequency counts. The 94 GHz 
counts are derived by making the additional assumption that 
all flat-spectrum sources have a simple power-law spectrum 
with a = -0.1 above 20 GHz (De Zotti et al. 2005). 

The De Zotti et al. (2005) model counts at 20 and 
30 GHz are calculated over a wide range in flux density 
(0.01 mjy to 10 Jy), though De Zotti et al. (2005) note that 
there are large uncertainties below ~0.1 mjy. Both Sunyaev- 
Zeldovich signals and free-free emission from proto-galactic 
plasma are likely to contribute significantly at sub-mjy flux 
levels, but there are currently almost no observational data 
available to constrain the models. 

Waldram et al. (2007) have recently estimated the 
source counts at 90 GHz based on the 15 GHz source counts 
from their 9C survey (Waldram et al. 2003) combined with 
follow-up 22 and 43 GHz VLA observations of a flux-limited 
sample of 121 sources selected at 15 GHz. The Waldram et 
al. (2007) source-count predictions are valid for the flux- 
density range between 10 and 200 mjy, since their 15 GHz 
survey area contains very few strong sources. 

Our derived 95 GHz source counts are in good agree- 
ment with the counts published by Waldram et al. (2007) in 
the flux density range where both data sets are valid (80- 
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Figure 9. Differential source-count predictions from AT20G 
20 GHz counts scaled to 95 GHz (this paper); Dc Zotti ct al. 
(2005) model counts at 94 GHz; Waldram et al. (2007) at 94 GHz 
and Holdaway et al. (1994) at 90 GHz. The vertical solid line rep- 
resents the lower bound of validity for our counts (82 mjy) , and 
the vertical dashed line is the upper bound of validity for the 
Waldram et al. counts (200 mjy). 
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Figure 10. Differential 95 GHz source counts normalized to the 
Euclidean case for the same four estimates shown in Figure 9. 
The flux— density scale now extends down to 1 mjy. As before, the 
vertical solid line represents the lower bound of validity for our 
counts. The two vertical dashed lines show the lower and upper 
bounds of validity for the Waldram et al. (2007) predictions. 
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>0.1 


Jy 
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± 1132 


120 ± 64 


De Zotti ct al. (2005) 


94 


4457 




288 


Holdaway ct al. (1994) 


90 


4400 




178 



Table 13. Comparison of the all-sky integrated source counts 
above a limiting flux density Sii m at the frequency v between our 
predictions and the Dc Zotti ct al. (2005) model and Holdaway at 
et. (1994) estimates. Our quoted errors have been derived from a 
Monte Carlo simulation based on the 20 GHz source counts and 
20-95 GHz spectral-index distribution, and arc dominated by the 
small-number statistics in the spectral-index distribution. 



200 mjy; see Figure 10). This is encouraging, as the two 
derivations are based on completely independent data sets. 

In contrast, our 95 GHz counts are significantly lower 
than the values derived by Holdaway et al. (1994) for all 
but the brightest sources, and are also lower than the De 
Zotti values at all flux densities. In the case of the De Zotti 
counts, the discrepancy almost certainly arises because De 
Zotti assumed a power-law spectral index of —0.1 at 20- 
94 GHz, whereas our 95 GHz data show that the median 
spectral index of AT20G sources in this frequency range is 
significantly steeper (~ —0.4, see Table 9). 

For the Holdaway et al. (1994) counts, the situation is 
more complex. As they note in their paper, the exact form 
of the inverted-spectrum tail of their 8-90 GHz spectral- 
index distribution (which is poorly defined because of the 
small number of sources observed in this region) is very im- 
portant because it allows a small fraction of the numerous 
weak 5 GHz sources to be boosted up to very high flux den- 
sities at 90 GHz. Holdaway et al. (1994) do not discuss the 
effects of variability over the several-year interval between 
their 8 and 90 GHz datasets, but this could also scatter vari- 
able objects into the tail of their spectral-index distribution. 
Since our results and those of Waldram et al. (2007) agree 
well and are also based on higher-frequency (15-20 GHz) 
source-count data, we conclude that Holdaway et al. (1994) 
have over-estimated the surface density of 95 GHz sources 
by up to a factor of two. 



5.5 Integrated source counts at 95 GHz 

As a consistency check, we can use the differential source 
counts derived in the previous section to estimate the num- 
ber of extragalactic radio sources brighter than a certain 
Siim over the whole sky at 90-95 GHz: 



N(S m > S lin 



4"/T 



f 



n(S) dS 



(4) 



Our results for Sn m = 0.1 and 1.0 Jy are reported in Ta- 
ble 13 together with predictions from De Zotti et al. (2005) 
and Holdaway et al. (1994). Since the Waldram et al. (2007) 
counts are only valid for sources below 200 mjy, integrated 
source counts cannot be derived from their study. 

Our derived source numbers are significantly lower than 
those derived by both De Zotti et al. (2005) and Holdaway 
(1994). This is a reflection of the lower differential source 
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counts found in our study, as discussed in the previous sec- 
tion. The actual number of sources stronger than 1 Jy over 
the whole sky at 90-95 GHz is still poorly determined, but 
the WMAP catalogue (Hinshaw et al. 2007) contains 106 
sources which lie above this flux level at 94 GHz. The ATCA 
calibrator catalogue at declinations south of 0° contains 60 
extragalactic sources stronger than 1 Jy in the 86-100 GHz 
band, implying roughly 120 such sources over the whole sky. 
Our derived value in Table 13 therefore appears reasonably 
consistent, within the errors, with what is observed. 

5.6 Confusion noise in CMB surveys at 95 GHz 

The 95 GHz source counts can also be used to compute the 
confusion noise fconf an d the angular power spectrum 5Ti of 
the temperature fluctuations due to discrete point-like extra- 
galactic radio sources, seen as a small angular scale contam- 
ination to the Cosmic Microwave Background anisotropies. 
This study is complementary to the work described in the 
previous section, because once the brightest sources in the 
90-100 GHz sky have been identified and removed from the 
CMB maps, the remaining fainter sources act as a confusing 
background. The formalism to derive a 2 and 6T t is based on 
Cleary et al. (2005) and briefly described below. 

Once a limiting flux density (Su m ) for a survey at a cer- 
tain observing frequency has been determined, the contribu- 
tion to the temperature fluctuations from faint sources below 
that flux density limit has two components: one C£° lsson due 
to the Poisson-distributed sources and a second term due to 
clustering C s c ^ stcring : 

Tp2 . — f /^-ipoisson . ^clustering \ /r\ 

CMB L-src — (, ) (,O src + O src ) (0) 

Csrc is the radio source angular power spectrum and 
dB v _ 2fc k Tqmb ,2 x A exp x 

~dT ~ ~& ( h ' (exp*-l) 2 [ ' 

is the conversion factor between fluctuations in the 
background intensity and in antenna temperature, k here 
is the Boltzmann's constant and x = hv/kTcMB, where 
Tcmb = 2.726 K is the CMB mean temperature. The term 
C s p ° isson given by the formula: 

ClT° n = a c 2 onf = /" Shm S 2 n(S) dS (7) 
Jo 

is also known as confusion noise for Poisson-distributed 
sources. This quantity is easy to estimate starting from ex- 
perimental differential source counts n(S). 

The term for source clustering (Scott & White 1999) in 
equation (5) is given by: 

/"Slim 

^clustering = WJ ( / g n(S .) dS) 2 (g) 

Jo 

where u>i is the Legendre transform of the angular two- 
point correlation function of radio sources uj(9). Lacking any 
clear estimate of the clustering properties of extragalactic 
radio sources at the frequency of 90 GHz, we assume for 
now that this term is negligible. 

In order to compare our prediction with others we con- 
ventionally expressed the point-source angular power spec- 
trum in terms of the Legendre multiple coefficient /: 



Author 


V 


A 




(GHz) 


(UK) 


This paper 


95 


0.78 


De Zotti ct al. (2005) 


94 


1.12 


Toffolatti et al. 1998 


100 


1.0 


Giommi et al. 2006 


94 


1.3 



Table 14. Comparison between the normalization coefficient A of 
temperature fluctuation angular power spectrum due to the back- 
ground of extragalactic discrete point-like sources fainter then 1 
Jy at the observing frequency v. 



ArT1 2 rji2 1(1 + 1) Csrc , I ^2 /„N 

AT src = T CMB — « const, x ( — ) (9) 

so that the STi is simply given by: 

ST, = (ATL) 1/2 ^xl (10) 

The normalization coefficient A of the angular power 
spectrum 5Ti obtained from our 95 GHz predicted counts 
is shown in Table 14 together with the model by Toffolatti 
et al. (1998), the model by De Zotti et al. (2005) which 
includes the contribution of flat-spectrum radio quasars, BL 
Lacs and steep-spectrum radio galaxies and the estimate 
by Giommi et al. (2006) for the blazar population only. In 
our computation of STi we assumed that the 95 GHz counts 
keep the same functional form for 695 < 82 mjy, which is 
the effective 95 GHz limit for our target sample. 



6 SUMMARY AND FUTURE WORK 

In this study, we have shown that the ATCA can measure 
robust 95 GHz (3 mm) flux densities for relatively faint con- 
tinuum sources in a few minutes of integration time, pro- 
vided that careful attention is given to calibrating the tele- 
scope pointing. Our ATCA measurements at 20 and 95 GHz 
appear to be on the same flux scale as the WMAP cata- 
logue (Hinshaw et al. 2007) to within the errors. The ar- 
rival in 2008 of the new ATCA broadband correlator, which 
has a 2 GHz bandwidth compared to the current value of 
128 MHz, will significantly increase the continuum sensitiv- 
ity at 95 GHz and should make it possible to measure 95 GHz 
sources down to ~ 10 mjy with integration times similar to 
those used here. 

We have made a new estimate of the 95 GHz radio 
source counts for sources stronger than 80mJy, and show 
that some previous studies have over-estimated the number 
of these sources by up to a factor of two. In particular, the 
confusion level in the Planck beam at 95 GHz is likely to be 
significantly lower than previously estimated. 

Fewer than half the sources in our 95 GHz sample cur- 
rently have measured redshifts, and obtaining good-quality 
optical spectra for the whole sample would be valuable both 
to study their rest-frame radio spectra and to test whether 
the brightest 95 GHz sources are relativistically beamed in 
the optical, as suggested in §4.8. 

We plan to make new measurements of a complete sam- 
ple of sources from the AT20G Bright Source Sample (Mas- 
sardi et al. 2008) at 20, 40 and 95 GHz in late 2007, using 
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the new ATCA 7 mm receiver system. These observations, 
which will also include polarization measurements, should 
provide further insights into the properties of the strongest 
95 GHz sources. 
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APPENDIX A: NOTES ON INDIVIDUAL 
SOURCES 

0058-3234: The 6dFGS spectrum (Jones et al. 2004) shows 
a bright featureless continuum and no redshift could be mea- 
sured. 

0111-4749: The optical ID is a 2dFGRS (Colless et 
al. 2001) galaxy at z=0.154 whose spectrum shows strong, 
narrow emission lines. This is one of the few AT20G sources 
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in our sample which was undetected at 95 GHz. The ra- 
dio spectrum rises steeply from 0.8 to 20 GHz, and so must 
peak somewhere above 20 GHz. The 843 MHZ SUMSS im- 
ages shows extended emission with a deconvolved largest 
angular size of 83 arcsec at PA 36° . 

0136—4044: This appears to be the same object as 
PMNJ 0136-4044, though the catalogued PMN position is 
significantly offset from the more accurate AT20G radio po- 
sition. There is a faint (stellar) optical ID at the AT20G 
position. 

0142—4206: This source appears extended in the 
843 MHz SUMSS image (deconvolved major axis 59 arcsec 
at PA 25°). 

0245—4459: The 6dFGS spectrum shows strong, broad 
Balmer emission lines typical of a Seyfert 1 galaxy or QSO. 

0320—3837: This source is so close to a bright star 
that no optical is ID possible. 

0428-4357: Blank field on the blue DSS image. No 
known optical ID. 

0433—4502: One of the few sources in our sample 
which shows an inverted radio spectrum at all frequencies 
from 0.8 to 95 GHz. The DSS optical ID is a 19th magnitude 
galaxy whose redshift is currently unknown. 

0440—4731: This source is well below the limit of 
the SUMSS catalogue at 843 MHz, but examination of the 
original SUMSS image shows a marginal (3 a) detection of 
a 2.5 mjy source at the AT20G position. There is no optical 
ID, and this is one of the few AT20G sources which does not 
appear in any previous optical or radio catalogue. 

0523—2955: This source is extended (or perhaps 
confused) in the SUMSS image (deconvolved major axis 
47 arcsec at PA 141°). The NVSS image also shows extended 
emission with largest angular size of 56 arcsec. 

0608—3041: This source is clearly extended north- 
south in the SUMSS image (deconvolved major axis 
92 arcsec at PA 1°), and an extended structure is also clearly 
seen in the AT20G 5 GHz image even though the image qual- 
ity is poor in this snapshot observation. The 1.4 GHz NVSS 
resolves the source into two components separated by about 
70arcse. At 20 GHz the source is unresolved. The 6dFGS 
spectrum shows strong, broad Balmer emission lines typi- 
cal of a Seyfert 1 galaxy or QSO. At z=0.2237, this is a 
relatively nearby object. 

1117-4837: This source is extended in the SUMSS 
image (deconvolved major axis 47 arcsec at PA 75°). 

1427—3305: The 6dFGS spectrum shows a featureless 
continuum, and no redshift could be measured. The optical 
counterpart is stellar on the blue DSS image, but is fuzzy 
(and classified as a galaxy) on the red image, suggesting that 
this may be a relatively nearby BL Lac object. 

1454—3747: The 6dFGS spectrum shows strong, 
broad Balmer emission lines typical of a Seyfert 1 galaxy 
or QSO. 

2031-3647: The 6dFGS spectrum shows emission 
lines of CIII] and Mgll. 

2113-3838: The 6dFGS spectrum shows a strong blue 
continuum with a single broad emission line of Mgll. 

2113-3838: The 6dFGS spectrum shows strong, 
broad emission lines of H7 and Mgll. 



Table 4. Table of 2005 observations. 'C after the source name in the first column indicates an ATCA calibrator, * indicates a source 
within 10° of the Galactic Plane (optical B magnitudes are not listed for these objects). 



2005 2005 2004 2004 2004 NVSS SUMSS 

Name Radio position Alt. name Bj T S95 S20 S20 ± S8.6 ± S4.8 ± S1.4 ± S843 ± z ID 

(J2000) mag Dayl Day2 mjy mjy mjy mjy mjy mjy 



(a) Sources with inverted radio spectra at 5-20 GHz 
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9.9 






0532-3848 


05 32 


02.14 


-38 48 


54.3 


PMNJ 0532-3848 


19.40 


2 


197 


192 


271 


235 


10 


190 


3 


131 


2 


168.6 


5.1 


139.8 


4.4 






0538-4405C 


05 38 


50.35 


-44 05 


08.6 


WMAP 148 


15.77 


2 


8051 


8076 


7328 


5294 


224 


4231 


60 


3800 


35 






3524.0 


105.8 


0.894 


Q 


0552-4522 


05 52 


10.00 


-45 22 


24.6 


PMNJ 0552-4522 


18.36 


2 


117 


133 


227 


237 


10 


173 


3 


143 


2 






248.3 


7.7 






0559-4529 


05 59 


11.53 


-45 29 


40.5 


PKS 0557-454 


18.34 


2 


546 


576 


797 


606 


27 


422 


6 


321 


3 






407.9 


12.3 


0.687 


Q 


0608-3041 


06 08 


40.59 


-30 41 


33.8 


PKS 0606-306 


17.83 


2 


82 


81 


124 


99 


3 


73 


5 


67 


5 


257.5 


8.7 


325.5 


12.2 


0.2237 


Q 


1102-4404 


11 02 


04.88 


-44 04 


22.3 


PKS 1059-438 


20.13 


2 


463 


493 


488 


765 


20 


703 


19 


566 


11 






111.0 


3.6 






1109-4815 


11 09 


18.96 


-48 15 


18.7 


PMNJ 1109-4815 


21.29 


1 


315 


335 


342 


335 


9 


157 


7 


170 


3 






193.5 


5.9 






1117-4838 


11 17 


20.01 


-48 38 


09.8 


PMNJ 1117-4838 


21.28 


2 


114 


135 


146 


204 


5 


127 


4 


96 


2 






97.1 


4.7 






1149-3723 


11 49 


01.83 


-37 23 


44.0 


PMNJ 1149-3724 


19.78 


2 


<90 


<73 


50 


73 


3 


77 


3 


73 


2 


71.8 


2.2 


82.0 


2.7 






1427-3305 


14 27 


41.55 


-33 05 


32.0 


WMAP 193 


20.46 


2 


1471 


1471 


871 


770 


21 


606 


17 


363 


10 


161.0 


4.9 


187.1 


6.2 






1454-3747C 


14 54 


27.57 


-37 47 


32.8 


PKS 1451-375 


16.36 


2 


947 


993 


1151 


1891 


49 


792 


22 


656 


16 


942.7 


28.3 


2045.0 


61.4 


0.3142 


Q 


1616-4018* 


16 16 


59.63 


-40 18 


04.3 


PMNJ 1617-4017 






95 


86 


120 


100 


3 


72 


2 


66 


2 






48.7 


1.9 






1624-4206* 


16 24 


58.72 


-42 06 


57.1 


PMNJ 1625-4207 






62 


56 


68 


75 


3 


48 


2 


35 


1 






33.4 


1.7 






1710-4718* 


17 10 


43.99 


-47 18 


20.3 








<93 


<70 


62 


76 


4 


28 


1 


16 


1 














2113-3838 


21 13 


53.91 


-38 38 


27.3 




18.45 


2 


<98 


63 


94 


101 


4 


71 


5 


59 


2 


44.6 


1.5 


57.2 


2.0 


0.850 


Q 


2151-3027C 


21 51 


55.56 


-30 27 


53.3 


PKS 2149-306 


17.95 


2 


621 


688 


2108 


1846 


58 


1314 


73 


1302 


36 


1243.6 


37.3 


1088.0 


32.7 


2.345 


Q 


2259-4818 


22 59 


26.26 


-48 18 


25.3 


PMNJ 2259-4818 


18.35 


1 


60 


82 


79 


72 


3 


69 


17 


47 


2 






54.5 


2.2 


0.1708 


G 



Table 4. Table of 2005 observations (continued). 











2005 


2005 


2004 




2004 




2004 




NVSS 




SUMSS 








Name 


Position 


Alt. name 


S95 


S20 


S20 


± 


S8.6 


± 


S4.8 


± 


SI. 4 


± 


S843 


± 


2 


ID 




(J2000) 




Dayl 


Day2 


mjy 


mjy 




mjy 




mjy 




mjy 




mjy 








(b) Other sources observed but later removed from the 'inverted spect 


rum' sample 


























0459-4528 


04 59 21.11 


-45 28 45.7 


PMN J 0459-4529 


<67 


<78 


45 


76 


2 


60 


2 


70 


2 






244.3 


7.5 






0537-3357 


05 37 32.14 


-33 57 55.2 




<86 


<58 


61 


74 


3 


50 


1 


54 


1 


32.3 


1.1 


58.8 


2.8 






0643-3314 


06 43 34.55 


-33 14 30.4 


PKS 0641-331 


238 


239 


212 


181 


5 


68 


5 


87 


4 


220.8 


6.6 


296.1 


9.2 






1109-3732 


11 09 57.71 


-37 32 20.0 


NGC 3557 


<118 


<88 


83 


52 


4 


23 


2 


77 


2 


366.4 


11.0 


567.4 


17.2 


0.0102 


G 


1132-3227 


11 32 27.53 


-32 27 14.1 




65 


68 


97 


85 


3 


106 


3 


87 


2 


65.2 


2.0 


67.5 


2.8 






1711-4552* 


17 11 05.99 


-45 52 54.1 




<115 


46 


93 


72 


4 


67 


2 


91 


2 






65.5 


2.6 






1742-3435* 


17 42 20.69 


-34 35 15.2 




81 


95 


144 


136 


7 










21.9 


0.9 


22.2 


1.7 






1745-4605* 


17 45 35.32 


-46 05 24.0 


IC 1266 


<101 


<83 


102 


68 


4 


56 


2 


65 


2 






113.9 


4.0 


0.000 


PN 


1913-3630 


19 13 21.03 


-36 30 19.7 


PMNJ 1913-3630 


98 


104 


119 


166 


13 










105.8 


3.2 


119.8 


43.8 






1937-3958C 


19 37 16.22 


-39 58 01.6 


PKS 1933-40 


1199 


1291 


1669 


1759 


45 










953.4 


28.6 


937.1 


28.2 


0.975 





1956-3225C 


19 56 59.57 


-32 25 45.1 


PKS 1953-325 


550 


643 


1000 


709 


51 










473.6 


14.2 


474.1 


14.4 


1.242 


o 


2000-4749 


20 00 12.98 


-47 49 51.7 




<83 


<85 


73 


79 


4 














7.4 


1.0 






2003-4336 


20 03 13.24 


-43 36 33.1 


PMNJ 2003-4336 


<131 


55 


96 


77 


4 














59.8 


2.0 






2025-3931 


20 25 42.54 


-39 31 40.3 


PMNJ 2025-3931 


<118 


<100 


52 


68 


3 










52.2 


1.6 


67.6 


2.3 






2031-3647 


20 31 49.14 


-36 47 42.1 




<108 


<75 


70 


62 


3 










33.2 


1.1 


28.6 


1.5 


1.5683 


o 


2039-5001 


20 39 48.79 


-50 01 12.1 


MRC 2036-501 


92 


89 


102 


105 


5 














524.6 


19.7 




G 


2056-4714C 


20 56 16.42 


-47 14 48.1 


WMAP 208 


1143 


1086 


1993 


1711 


55 














2138.0 


64.2 


1.489 


O 


2223-3137 


22 23 21.69 


-31 37 03.0 


PKS 2220-318 


485 


476 


642 






447 


10 


329 


6 


220.7 


6.6 


274.7 


8.4 






2235-4835C 


22 35 13.27 


-48 35 58.8 


WMAP 206 


2605 






1994 


75 


2065 


518 


1207 


35 






1216.0 


36.5 


0.510 





2235-3629 


22 35 54.86 


-36 29 02.1 




268 


272 


246 






156 


39 


120 


4 


68.2 


2.1 


74.4 


2.7 






2245-4931 


22 45 00.24 


-49 31 48.2 


PKS 2242-498 


<114 


<81 


84 


97 


4 


51 


23 


88 


3 






310.3 


9.4 






2248-3235C 


22 48 38.81 


-32 35 50.8 


PKS 2245-328 


1166 


1170 


1655 






826 


207 


410 


12 


708.0 


21.2 


710.5 


21.4 


2.268 


Q 


(c) Strong sources from the ATCA calibrator list 
































0004-4736C 


00 04 35.65 


-47 36 19.2 


PKS 0002-478 


893 






868 


16 


970 


59 


900 


30 






932.9 


28.0 






0012-3954C 


00 12 59.89 


-39 54 26.4 


WMAP 202 


734 


716 


1420 


1609 


33 


2009 


121 


2015 


67 


465.0 


14.0 


839.9 


25.3 






0051-4226C 


00 51 09.50 


-42 26 32.6 


PKS 0048-427 


708 






1462 


32 


1544 


93 


1267 


42 






689.8 


20.8 


1.749 


Q 


0440-4333C 


04 40 17.16 


-43 33 08.4 


WMAP 147 


460 


439 


1533 


1947 


51 


2961 


74 


3681 


28 






5858.0 


175.8 


2.863 


Q 


0631-4154C 


06 31 12.00 


-41 54 27.1 


PKS 0629-418 


165 


162 


399 


517 


17 


623 


35 


708 


40 






684.9 


20.6 


1.416 


Q 


0648-3044C 


06 48 14.17 


-30 44 19.2 


PKS 0646-306 


780 


851 


929 


744 


24 


713 


40 


772 


43 


1295.2 


1.0 


985.3 


29.7 


0.455 


Q 


1107-4449C 


11 07 08.58 


-44 49 06.8 


WMAP 166 


645 






1674 


40 


2174 


59 


2703 


50 






2516.0 


75.5 


1.598 


Q 


1454-4012C 


14 54 33.19 


-40 12 31.5 


PKS 1451-400 


222 


246 


462 


526 


15 


571 


16 


719 


16 


960.9 


28.8 


1210.0 


36.3 


1.810 


Q 


2158-3013C 


21 58 51.93 


-30 13 31.0 


PKS 2155-304 


310 


327 


388 


350 


10 


377 


21 


369 


10 


424.2 


12.7 


479.2 


14.6 


0.116 


Q 



Table 5. Tabic of 2006 observations. 



Name 






Radio position 
(J2000) 


Alt. name 




T 


2006 
S95 
mjy 


2006 
S20 
mjy 


2004 
S20 
mjy 


± 


2004 
S8.6 
mjy 


± 


2004 
S4.8 
mjy 


± 


NVSS 
S1.4 
mjy 


± 


SUMSS 
S843 
mjy 


± 


z 


ID 


oooft 
uuuo— 


494^ 
-4Z40 


nn 


06 


19.72 


-42 45 18.3 


pt<^ nnn^ A^n 

rj\c UUUO — ^±OU 


1 8 97 


9 
Z 


oo 


1 48 


1 S3 
loo 


A 
1 


1 8Q 

ioy 


1 1 
1 1 


1 QQ 

iyy 


7 






94ft 
Z4D.U 


7 7 


1 780 
1 . ( oU 




001 9 
uuiz— 




uu 


12 


59.89 


-39 54 26.4 


WM A P 9n9 
VV lvl/\r ZUZ 


1 8 33 
lo.oo 


9 
Z 


4Q£ 

^yo 


Q&4 

yo4 


1 ROQ 

louy 


33 
OO 


900Q 

zuuy 


191 
1Z1 


901 ^ 
ZU10 


(17 


4Q3 ft 

^tyo.o 


17 4 
1 ( .4 


83Q Q 

ooy.y 


9^ 3 
ZO.O 






0094 


-4202C 


nn 
uu 


24 


42.95 


-42 02 03.5 


PT<T^ nn99 493 
1 l\ i J uuzz — ^zo 


99 7R 


2 


'i-O 


9fi7 

zu / 


987 
ZO I 


7 


QQ9 

yyz 


fiO 
UU 


1 Q^fi 

iyou 


OO 






1 770 
1 / I u.u 


^3 1 
OO. 1 


Q37 

u.yo i 


n 


0030 

UUOU — 


4994 


nn 
uu 


30 


17.51 


-42 24 46.3 


PT^t; nn97 Ai(\ 

rj\o UUZ l — ^zo 


lO IK 

lo.lO 


9 

z 


974 


3SQ 

ooy 


47S 


1 
1U 


^1 1 

Oil 


31 
01 


M 1 

Oil 


1 7 
1 I 






494 8 
4Z4.0 


19 8 
1Z.O 


4Q^ 

u.4yo 


n 


UU4U 


3943 
-OZ40 


nn 
uu 


40 


17.63 


-32 43 27.6 


PMNinnzin 39zl^ 

.r 1V11NJ UU*±U — OZ^±0 


91 

zi .oy 


9 


128 


98^ 
ZOO 


31 4 

Oil 


r: 
o 


97fi 
z / u 


1 7 
1 t 


231 


Q 

o 


91 7 
Z1U. i 


ft 3 

u.o 


1 78 9 
1 ( o.z 


^ ft 
o.u 






00/10 
UU4U — 


399^ 

-ozzo 


nn 
uu 


40 


30.73 


-32 25 19.9 


PT^t; nn^R 39^ 
r rvo UUoo — OZD 


91 41 
Z± .*±1 


9 
Z 


oy 


1 «7 
ID I 


1 Q3 

iyo 


A 
*± 


900 
ZUU 


1 9 
1Z 


1 Q4 

ly^ 


7 


1114 
1 1 1.4 


3 4 
0.4 


1 38 1 
lOo. 1 


4 ^ 
4.0 






00^8 

UUOO 




nn 
uu 


58 


02.22 


-32 34 20.5 


PT^nn^ *398 

-T -Tvo UUOO OZO 


90 03 

zu.uo 


2 


140 


1 Qfi 

iyo 


397 
OZ / 


fi 
U 


30Q 

ouy 


1 Q 

iy 


9^ 
ZUO 


1 
1U 


18^0 
loo.y 


^ ft 
O.O 


184 1 


^ 8 
O.O 






00^8 


3347 
004 t 


nn 
uu 


58 


15.70 


-33 47 57.1 


.r IVliNJ UUOO — oo*±l 


90 ^4 


9 
Z 


1 03 
1UO 


1 3Q 

ioy 


1 74 
1 { ^± 


A 
*± 


91 7 
Zl ( 


1 3 
lO 


914 
Z14 


7 


1 9Q 9 

izy.z 


3 Q 

o.y 


118 4 
1 lo.4 


3 8 
O.o 






01 Oft 


4034^ 
"4U04V_< 


m 

u ± 


06 


45.11 


-40 34 19.5 


WM Ap 171 


1 R fiO 


2 


91 QQ 
z±yy 


41 fiO 


2146 


7Q 

/ y 


1 ^40 

lo^y 


38 
Oo 


1 OQQ 

iuy y 


77 






ft^O Q 

uou.y 


1 Q ft 

iy.u 


^84 
U. Oo4 


n 
m 


01 1 8 


41 09 
-41UZ 


m 
u 1 


18 


45.27 


-41 02 30.6 


PT^t; ni i £ A^ i 


1 Q 44 


9 
Z 


£9 
oZ 


loD 


1 79 
1 1 z 


D 


91 1 
Z 1 1 





998 
ZZo 


1 Q 

iy 






1 090 4 
1UZU.4 


30 7 
OU. I 






01 33 
U 1 OO — 


4430 
-440U 


ni 
u 1 


33 


00.87 


-44 30 43.3 


PT^t; ni ^n AA7 
rj\o UloU — <±*± i 


J> ZZ 


n 
U 


SI 


90^ 
ZUO 


91 ^ 
Z10 


Q 
O 


90.^ 

zyo 


Q 

o 


499 

^iZZ 


31 
Ol 






4^^ 1 
400. 1 


13 8 
lo.o 






01 34 
U 1 04 — 


■JO/IIP 


ni 
u 1 


34 


32.14 


-38 43 33.7 


ivirxo uioz — ooy 


1 7 7fi 


9 

z 


D^O 


SOQ 

ouy 


DoO 


93 
ZO 


^07 

ou / 


1 9 
lz 


438 
^tOo 


31 
Ol 


^ftQ 1 

ouy. i 


17 1 
11.1 


ft^Q Q 

ooy .y 


1 Q 8 


9 1 40 
Z. 14U 


rv 


01 3fi 
U i. oo 


-4044 


u ± 


36 


02.88 


-40 44 50.4 


PA/TN T ni 3fi zLnzlzl 


91 3fi 
Zl .ou 


9 


1 QO 

iyu 


974 

Z / 4 


3^0 
OOU 


11 


973 
z / o 


7 


226 


1 ft 
1U 






115 1 


3 ft 
O.O 


ft4Q 
u.U4y 


o 


01 44 
U144 — 


3Q38 


ni 
u 1 


44 


54.10 


-39 38 10.8 


PMW T m ztzl 3Q37 


91 


9 
Z 


7S 
I o 


1 33 
lOO 


900 
ZUU 


7 


998 
zzo 


a 
D 


99Q 

zzy 


1 ft 
lO 


8Q 4 


9 7 
Z. i 


89 9 
oz.z 


9 7 
Z. i 






01 ^ 


4048P 


ni 

u ± 


55 


37.07 


-40 48 42.4 


pt<t^ m zli n 

I l\ i j u±oo — ftlU 


1 Q 31 




125 


40^ 


426 


14 


800 

ouu 


90 
zu 


1252 


88 
oo 






2426 


79 8 
( z.o 


226 


r; 

VT 


01 ^7 

uio / — 


389Q 

-oozy 


ni 

U 1 


57 


19.72 


-38 29 44.4 




90 41 


9 
Z 


oD 


1 ^4 


939 
ZOZ 


Q 
o 


1 

ioy 


A 
1 


Q1 

y i 


7 


9Q 9 

zy.z 


1 
1 .u 


97 ft 

z t .o 


1 7 
1. ( 






090ft 
UZUO — 


3094 
-OUZ4 


n9 


06 


43.28 


-30 24 58.4 


PT^t; n9n4 ^nfi 

1 l\ o UZU^t — OUO 


99 fi9 
ZZ.DZ 


9 
Z 


00 


1 (^1 
1D1 


1 80 
loU 


D 


349 

o^tz 


A 
1 


4QQ 

^yy 


7 


844 3 

0*4:4.0 


9^ 3 

zo.o 


Q87 1 

yo i . i 


9Q 7 

zy. ( 




P 
V T 


091 ft 
UZIO — 


3947 
-OZ4 i 


ft 9 
UZ 


16 


48.19 


-32 47 40.6 


PT^t;n9i4 ^n 

1 l\ O UZl^t — oou 


lo.oo 


9 
Z 


970 
Z / U 


4^^ 

^too 


^40 
04U 


1 ^ 
lO 


^04 
OU4 


O 


330 

oou 


O 


1 ft 1 ^ 9 
1DO.Z 


^ 

o.u 


KO 1 
10Z. 1 


4 Q 

4.y 


1 331 
l.OOl 


o 


0997 

uzz / — 


309ft 
OUZO 


no 


27 


40.55 


-30 26 04.2 


PT^t;n99^ ^nfi 
1 i\o uzzo — OUO 


1 S 70 

lo. / y 


9 
Z 


30S 
OUO 


904 


4^9 
^±OZ 


1 3 
lO 


OOO 


1 


31 O 
01U 


'1 


9Q3 ft 

zyo.u 


8 8 
O.O 


9^1 Q 

zoi .y 


8 1 

O. 1 


303 
U.OUO 


o 
H 


0931 


474ft 
-4 i 40 


ft 9 

uz 


31 


11.77 


-47 46 12.0 


PMN T n931 Zl7z1fi 


1 7 87 


9 
Z 


41 R 


S34 


77^ 
I i 


1 

iy 


OOO 


a 
D 


3(^3 
OOO 









49 8 
4Z.O 


1 ft 
1 .o 






0938 


3031 
OUOi. 


f)9 
uz 


38 


55.31 


-30 32 02.6 


PA/TNT n9' : t£ 3n^9 


1 R 90 


2 


1 Qfi 

iyu 


220 


1 73 
1 f o 


O 


1 8^ 
lOO 


2 


1 8fi 

lOU 


Q 
O 


1 ft 1 ^ 4 
IUO.4 


^ Q 

o.y 


910^ 
Z1U.O 


ft ^ 

O.O 






094^ 
UZ4D — 


44^0 

-44oy 


n9 


45 


54.07 


-44 59 39.5 


PT^t; noAA zL^9 
r 1 \ o uz^i^t — ^±OZ 


1 7 


i 

i 


9S7 
Zo i 


490 


^81 
Ool 


1 7 
1 / 


oyi 


D 


R88 
Ooo 


Q 

y 






1 Q79 

iy / z.u 


^Q 9 

oy .z 


980 
U.ZoU 




094ft 
UZ40 — 


4ft^1 
-4001 


f!9 
UZ 


46 


00.09 


-46 51 16.8 




1 £ 40 


9 
Z 


91 R 
Zlo 


3^3 
000 


3Q7 

oy / 


1 9 
1Z 


^i « 

OlD 


O 


fi77 
( I 


Q 

y 






1 ^90 
10ZU.U 


4^ ft 
40. 






0300- 


-4959 


03 


00 


27.80 


-49 59 11.8 




> 22 





<108 


141 


170 


6 


214 


8 


1/1 


3 






300.7 


9.1 






0317- 


-4414 


03 


17 


57.66 


-44 14 16.9 


ESO 248-G006 


15.74 


1 


112 


377 


258 


9 


449 


17 


611 


8 






1915.0 


57.5 


0.0761 


G 


0320- 


-3837 


0.3 


20 


46.43 


-38 37 28.4 


PKS 0318-388 






138 


162 


161 


5 


145 


6 


135 


2 


175.2 


5.3 


205.5 


6.3 






0330- 


-4014 


0.3 


30 


51.11 


-40 14 16.3 


PKS 0329-404 


22.40 


2 


186 


326 


311 


10 


470 


18 


577 


8 






938.4 


28.2 






0332- 


-4455 


0.3 


32 


44.10 


-44 55 57.1 


PKS 0330-450 


17.69 


2 


66 


195 


180 


6 


227 


9 


159 


2 






111.4 


3.5 


2.60 


Q 


0334- 


-4008C 


03 


34 


13.62 


-40 08 25.3 


WMAP 146 


17.51 


2 


353 


552 


1274 


42 


1124 


43 


943 


13 






986.5 


29.7 




Q 


0334- 


-3725 


0.3 


34 


15.49 


-37 25 43.2 


PMNJ 0334-3725 


17.08 


2 


322 


358 


293 


10 


319 


12 


313 


4 


221.2 


7.8 


171.4 


5.3 






0336- 


-3616 


03 


36 


54.11 


-36 16 06.0 


PKS 0335-364 


18.40 


2 


640 


805 


738 


22 


610 


23 


496 


7 


500.9 


15.0 


526.2 


15.9 


1.541 


Q 


0342- 


-3703 


0.3 


42 


05.41 


-37 03 21.0 


PKS 0340-37 


17.85 


2 


59 


236 


217 


7 


324 


13 


696 


9 


2059.6 


72.6 


3126.0 


93.8 


0.284 


Q 


0402- 


-3147C 


04 


02 


21.31 


-31 47 25.7 


PKS 0400-319 


20.33 


2 


308 


536 


467 


9 


587 


15 


632 


5 


681.3 


20.4 


506.5 


15.3 


1.288 


Q 


0407- 


-3303C 


04 


07 


33.92 


-33 03 45.3 


PKS 0405-331 


19.33 


2 


321 


692 


563 


14 


659 


16 


634 


5 


634.9 


19.1 


697.7 


21.0 


2.562 


Q 


0415- 


-4637 


04 


15 


58.69 


-46 37 49.2 


PKS 0414-467 


21.75 


1 


80 


198 


171 


5 


287 


7 


358 


3 






477.9 


14.4 







Table 5. Table of 2006 observations (continued). 



2006 

Name Radio position Alt. name Bj T S95 

(J2000) mjy 



2006 2004 2004 2004 NVSS SUMSS 

S20 S20 ± S8.6 ± S4.8 ± S1.4 ± S843 ± z ID 

mjy mjy mjy mjy mjy mjy 



0421- 


-4818 


04 


21 


27.49 


-48 18 05.3 


PKS 0419-484 


16.85 


2 


212 


321 


452 


13 


356 


9 


232 


2 






284.8 


8.6 


0.527 


O 


0422- 


-3315 


04 


22 


54.29 


—33 15 08.7 


PKS 0421—333 


21.72 


2 


93 


171 


179 


5 


286 


7 


427 


3 


950.5 


28.5 


1347.0 


40.5 






0439- 


-4522C 


04 


39 


00.83 


—45 22 22.6 


PKS 0437-454 


20.45 


2 


266 


442 


700 


17 


890 


22 


928 


7 






950.9 


28.6 






0441- 


-3340C 


04 


41 


33.67 


—33 40 03.7 


PKS 0439—337 


22.48 


1 


<107 


238 


269 


7 


541 


14 


799 


6 


1155.9 


34.7 


896.0 


27.0 






0441- 


-3300 


04 


41 


34.97 


—33 00 08.1 


PKS 0439-331 


19.86 


2 


83 


167 


160 


4 


175 


4 


186 


2 


140.9 


4.3 


233.0 


7.3 






0449- 


-4809 


04 


49 


39.59 


—48 09 40.9 


PMN J 0449-4809 


18.54 


2 


80 


284 


209 


6 


261 


7 


221 


2 






99.9 


3.1 






2003- 


-3251C 


20 


03 


24.04 


-32 51 42.5 


PKS 2000-330 


18.83 


2 


61 


401 


504 


22 










446.0 


15.7 


300.8 


9.2 


3.773 





2006- 


-3451 


20 


06 


46.31 


—34 51 14.6 


PMN J 2006-3450 


20.53 


2 


86 


169 


177 


8 










96.0 


2.9 


108.8 


3.5 






2009- 


-4849C 


20 


09 


25.45 


-48 49 53.9 


PKS 2005-489 


15.75 


1 


521 


737 


877 


42 














1299.0 


39.0 


0.071 


G 


2012- 


-4721 


20 


12 


09.52 


-47 21 20.2 


PMNJ 2012-4721 


20.74 


1 


52 


138 


175 


8 














298.9 


9.0 






2024- 


-3253 


20 


24 


35.69 


-32 53 36.8 


PKS 2021-330 


18.73 


2 


272 


438 


479 


23 










886.5 


26.6 


1201.0 


36.1 


1.471 





2027- 


-4726 


20 


27 


50.74 


-47 26 20.3 


PMNJ 2027-4726 


21.75 


1 


160 


254 


305 


15 














230.5 


7.0 






2056- 


-4714C 


20 


56 


16.39 


-47 14 48.1 


WMAP 208 


18.29 


2 


2453 


3246 


1711 


55 














2138.0 


64.2 


1.489 





2105- 


-4848 


21 


05 


01.20 


-48 48 46.8 


PKS 2101— 49 


17.31 


2 


173 


516 


468 


14 


687 


39 


735 


20 






1516.0 


45.5 






2108- 


-4846C 


21 


08 


50.80 


-48 46 21.8 


PKS 2105-48 


20.03 


2 


78 


212 


225 


7 


420 


24 


586 


16 






1459.0 


43.8 






2123- 


-3046 


21 


23 


48.69 


-30 46 06.2 


PMNJ 2123-3046 


19.18 


2 


57 


293 


179 


6 


271 


15 


300 


8 


147.2 


4.4 


155.8 


5.1 






2124- 


-4948 


21 


24 


29.12 


-49 48 06.8 


PMNJ 2124-4948 


19.46 


2 


<92 


74 


171 


6 


210 


12 


211 


6 






222.3 


6.8 






2145- 


-3142 


21 


45 


47.01 


-31 42 34.0 


PKS 2142-319 


20.82 


2 


129 


152 


165 


6 


184 


11 


169 


5 


144.3 


4.4 


188.1 


5.9 


0.66 





2151- 


-3027C 


21 


51 


55.58 


-30 27 53.8 


PKS 2149-306 


17.95 


2 


849 


1402 


1846 


58 


1314 


73 


1302 


36 


1243.2 


37.3 


1088.0 


32.7 


2.345 





2158- 


-3013C 


21 


58 


52.08 


-30 13 32.4 


PKS 2155-304 


12.35 


2 


903 


844 


350 


10 


377 


21 


369 


10 


489.3 


19.0 


479.2 


14.6 


0.116 


Q 


2229- 


-4051C 


22 


29 


18.61 


-40 51 32.5 


PKS 2226-41 


19.58 


2 


49 


277 


254 


10 


790 


198 


1128 


33 






4230.0 


129.6 


0.4462 


Q 


2230- 


-4416 


22 


30 


56.44 


-44 16 30.0 


PKS 2227-445 


19.46 


2 


350 


472 


472 


18 


628 


157 


423 


12 






520.1 


15.7 


1.326 


Q 


2235- 


-4835C 


22 


35 


13.29 


-48 35 58.7 


WMAP 206 


17.61 


2 


1076 


1775 


1994 


75 


2065 


518 


1207 


35 






1216.0 


36.5 


0.510 


Q 


2257- 


-3627C 


22 


57 


10.60 


-36 27 44.0 


IC 1459 


11.49 


1 


264 


621 






1208 


303 


1198 


35 


1279.7 


45.2 


966.2 


29.1 


0.0056 


G 


2300- 


-3536 


23 


00 


55.92 


-35 36 15.9 


PMNJ 2300-3536 


21.08 


2 


97 


174 


220 


3 


220 


9 


161 


3 


97.9 


3.0 


88.8 


3.1 






2302- 


-3718 


23 


02 


23.77 


-37 18 07.3 


PKS 2259-37 


> 22 





41 


240 


261 


4 


585 


24 


986 


13 


2719.8 


81.6 


4243.0 


127.3 






2302- 


-4533C 


23 


02 


40.04 


-45 33 57.7 


PKS 2259-458 


20.46 


1 


93 


148 


165 


3 


194 


4 


305 


4 






1059.3 


31.8 






2307- 


-3548 


23 


07 


37.27 


-35 48 28.7 




21.14 


1 


106 


118 


190 


3 


186 


8 


112 


2 


58.7 


2.2 


93.5 


3.1 






2308- 


-4557 


23 


08 


52.81 


-45 57 56.1 




19.39 


1 


199 


225 


286 


4 


212 


9 


163 


2 






214.3 


6.5 




G 


2313- 


-4729C 


23 


13 


51.91 


-47 29 11.7 


PKS 2311-477 


22.20 


1 


53 


181 


189 


3 


421 


18 


647 


9 






966.5 


29.0 






2314- 


-4455C 


23 


14 


09.38 


-44 55 48.9 


PKS 2311-452 


19.95 


2 


125 


439 


497 


7 


981 


41 


1391 


19 






2863.0 


85.9 


2.884 


Q 


2316- 


-4337C 


23 


16 


21.07 


-43 37 47.0 


PKS 2313-439 


19.82 


2 


124 


303 


319 


5 


523 


22 


682 


9 






1384.0 


41.6 


1.847 


Q 


2316- 


-4041 


23 


16 


46.90 


-40 41 20.7 


PKS 2314-409 


18.76 


2 


82 


247 


245 


4 


392 


16 


397 


5 






520.8 


15.8 


2.448 


Q 


2318- 


-4010 


23 


18 


06.88 


-40 10 06.5 


PKS 2315-404 


19.52 


2 


157 


231 


188 


3 


235 


10 


288 


4 






296.8 


9.1 


1.820 


Q 


2318- 


-4032 


23 


18 


09.29 


-40 32 06.7 


PKS 2315-408 


18.97 


2 


52 


166 


173 


3 


333 


14 


502 


7 






1329.0 


39.9 







